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1. Introduction 

The Bardeen-Cooper-Schrieffer (BCS) equation valid in the weak-coupling limit (Fermi- 
Fermi scattering length a — )■ —0) of a Fermi super-fluid (SF) has been fundamental for 
the study of low-temperature superconductors p. Later, a microscopic field-theoretic 
formulation of the Fermi SF has been proposed by Bogoliubov and de Gennes [1] . Eagles 
[2] suggested the possibility of extending the study of the Fermi SF to the domain 
of stronger coupling which initiated the study of the Fermi SF in the BCS to Bose- 
Einstein condensation (BEC) crossover |3] . As the strength of the attractive interaction 
between two spin-up and -down fermions in a (Cooper) pair forming the Fermi SF 
is increased (a < 0) the simple BCS SF turns into a complex Cooper-pair-induced 
strongly interacting SF at unitarity (a — >■ ±oo) and beyond (a > 0) these two fermions 
form a diatomic molecule (dimer) and undergo BEC. In this domain of interaction, the 
properties of the SF can largely be described by mean-field bosonic equations and in 
the so called BEC limit as a — )■ +0, the system is ideally described by the mean-field 
Gross-Pitaevskii (CP) equation [1] for dimers. Although the studies of the Fermi SF 
in the weak-coupling BCS (a — )■ —0) and BEC (a — )■ -|-0) limits as well as at unitarity 
(a — )■ ±oo) are well under control, its study in the full BCS-BEC crossover, as a varies 
from —0 to +0 through ±oo remains to be a challenging and formidable task |5]. 

After the first observation of a strongly interacting degenerate Fermi gas [6], 
recently, several experimental groups have observed and studied a trapped Fermi SF 
in the the BCS-BEC crossover composed of [7j and ^Li [HIE] atoms under controlled 
conditions. The crossover in the Fermi SF was attained by varying an external 
background magnetic field near a Feshbach resonance [10] which allows an experimental 
manipulation of the S-wave scattering length a. The study of the Fermi SF in the 
BCS-BEC crossover has gained new impetus [5l [TTl [121 [131 Ej after these experiments 
and the necessity of a proper theoretical formulation to study the trapped Fermi SF 
in the crossover cannot be overemphasized. Recently, we formulated fi5\ [T6] a density- 
functional (DF) equation [17] which produced results [TH |19] for energy of a spherically- 
symmetric trapped Fermi SF in the BCS-unitarity crossover in close agreement with 
Monte Carlo calculations [20]. Here we extend that DF equation to the whole BCS- 
BEC crossover and apply it to the study of stationary and non-stationary problems 
of a cigar-shaped [21] Fermi SF. Usually, the use of any microscopic formulation of a 
Fermi SF (including the Monte Carlo simulations) in a large system as encountered in 
laboratory is extremely complicated numerically and/or analytically and the advantages 
of a DF approach to Fermi systems have long been realized [17] . The DF approach with 
a small number of degrees of freedom is reasonably easily formulated and applied to 
large Fermi systems |17] . 

Before we detail the DF formulation for Fermi SF, we list the desired aspects 
and properties of the Fermi SF which are included in the present formulation, (i) A 
Fermi SF has been long described by a Galilean-invariant hydrodynamical formulation 
[22]. By showing complete equivalence of the present DF equation with a quantum 
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hydro dynamical equation, we preserve Galilean invariance of the present equation and 
include proper gradient (surface) correction [T5l [16] in the Thomas- Fermi approximation 
(or the local-density approximation) [5]. (Such surface correction was first introduced 
by von Weizsacker [23] in the study of properties of large nuclei.) (ii) There are model- 
independent results l25l [26] for energy per atom of a uniform Fermi SF near the BCS 
limit in terms of the scattering length a, as well as of a SF of dimers in the weak-coupling 
BEG limit, and these will be included in our formulation, (iii) A theoretical formulation 
of the the uniform Fermi SF at unitarity is particularly interesting as the only available 
scale of length at unitarity is n"^/'^, where n is the density, and all energies should have 
the universal form [271 EH] ~ h^n^^^ /m, where m is the mass of a pair. This aspect 
of the Fermi SF at unitarity is included in the present DF equation. (An interesting 
account of the universal behavior of the properties of a Fermi gas at unitarity has 
recently appeared [29].) There are other DF-type equations for a Fermi SF valid for the 
crossover [30l [3T| [32] and also at unitarity [11], but they fail on one or more accounts 
above. They all fail to satisfy criterion (i) above. The studies of [301 El] do not satisfy 
criterion (ii). It is not easy to realize if any of these properties are violated in the 
approximate solution of the Bogoliubov equation as reported in [33]. (The accurate 
solution of the Bogoliubov equation should not violate any of these properties.) 

In Sec. [2] we present and discuss the DF equation for Fermi SF in the BGS-BEG 
crossover. In Sec[3]we present the numerical results of the present formulation. In Sec l3.1l 
we discuss the stationary properties, such as size and linear density of a cigar-shaped 
Fermi SF, and compare our results with experiment [3l]. In Sec 13.21 we discuss the 
breathing oscillation of a cigar-shaped Fermi SF and compare the results for frequencies 
of radial and axial oscillation with experiment [35l [Ml [37]. Finally, in Sec. [4] we present 
some concluding remarks. 



2. DF equation for a Fermi SF in the BCS-BEC crossover 



We consider a Fermi SF composed of an equal number of fully paired spin-up and spin- 
down fermions. In the present study we start with the following DF equation [T8l [T9] 
which produced results for energy of a Fermi SF in a spherically-symmetric harmonic 
trap in close agreement with Monte-Garlo (MG) calculations [20] in the BGS-unitarity 
crossover [3] 



■V^ + f/(r) + yu(n,a) - i 



d 



^(r,t) = 0, 



(1) 



2m ' ' ^ ^ ■ r \ 1 J 

where U is the trapping potential, m is the mass of a pair (twice the atomic mass), 
n{n, a) is the bulk chemical potential defined by 



/i(n, a) 



Xix - X2X^ 



X < 0, 



(2) 



2Ef 1 - + f32X' 

where n = is the density of pairs, the gas parameter x = n^^^a, Ep = h'^(QTr'^ny^^ /m 
is the Fermi energy of the uniform gas. The normalization condition of the DF wave 
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function \l/ is / |\l/p(i^r = A^, where is the number of pairs. The Fermi- Fermi scattering 
length a is negative corresponding to attraction: —0 > a > — oo. The constants of ([2]) 
have the following values [H]: xi = 4.n2^'^ / {'in'^f'^ ,X2 = 2^/^300, f3i = and 
^2 = X2/(l ~ 0) where ^ is an universal constant (called Bertsch parameter) which 
determines the value of the bulk chemical potential at unitarity [271 l28l 138] : 

lim a(n, a) = 2Ep^. (3) 

a— ^±oo 

At unitarity the only length scale in a uniform Fermi SF is n~^/^, and from dimensional 
argument the chemical potential of the trapped SF fermions (and also the energy) have 
the above universal form [21]. The bulk chemical potential /i(n, a) of ([2]), has the 
following leading terms of dilute uniform gas in the BCS limit (a — —0), as obtained 
from the model-independent expression for energy per particle derived by Lee and Yang 
[MIES] 

u(n,a) ^ 47r2^/^ / i/q x , 

The value of the Bertsch parameter ^ of ([2D can be obtained from experiment on a 
trapped Fermi SF [5] or from the MC calculation of a uniform Fermi SF [271 [28] . There 
have been theoretical [271 [2S1 EB] and experimental [2] investigations which extracted 
the value of the constant ^ for a Fermi SF, and now there is a consensus [5] that its 
value should not depend on the specific Fermi atoms used. The experimental extraction 
of this constant from an analysis of trapped Fermi SF has yielded slightly different 
values ranging from 0.3 to 0.5 under different trapping conditions [H [32], whereas MC 
calculation in a uniform SF have yielded [27] the value 0.44. A recent average over 
experimental evaluations of the Bertsch parameter is ~ 0.38 [39] within about 2%. The 
values of this parameter are obtained from energy-entropy measurements of the ground 
state, sound speed, and cloud size. (The BCS mean-field [IQ] calculation yielded slightly 
different result: ^ = 0.59.) The spread of the experimental value of ^ under different 
trapping conditions seems to be well outside the range of experimental uncertainty. 
Hence a question naturally arises. Whether the value of ^ for a trapped Fermi SF 
should depend on trapping conditions and number of atoms. We explore this possibility 
in the present investigation. 

As the Fermi-Fermi attraction is increased beyond unitarity so that the scattering 
length a turns positive, the Fermi pairs, each composed of a spin-up and a spin-down 
fermion, are transformed into bound molecules (dimers) which undergo BEC. In the 
weak-coupling limit of weak interaction (BEC limit) the Fermi SF is described by the 
GP equation for dimers. A proper description of the BEC of dimers should be made 
in terms of the dimer-dimer scattering length a^. Using the zero-range approximation 
Petrov et al. [H] calculated the dimer-dimer scattering length in terms of the Fermi- 
Fermi scattering length a and obtained 

ad = 0.6a, (5) 

which we shall use in the present study. 
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In the BEC-unitarity crossover, as the dimer-dimer scattering length varies from 
+0 to +00, we shall adopt a recent beyond mean-field extension of the GP equation 
[12] for fundamental bosons valid in the weak-coupling limit to unitarity. For a BEG 
of dimers, the energy and chemical potential have the same universal form (|3]). (For 
fundamental bosons considered in |12], at unitarity, the constant ^ had a value different 
from ^ = 0.44, e.g. ^ = 0.73 extracted from (16) of [13].) We now generalize the bulk 
chemical potential fi{n, a) of ([2]) to the full BGS-BEG crossover. This can be achieved 
if we take n{n, a) of (II]) in the BEG-unitarity crossover (+00 > a,ad > +0) as 

t^{n,a) ^ 271 y+{l + v)ay^l'^ , . 

2Ef (67r2)2/3l + i,tt^3/2 + (l + ^)^^5/2' ^ ' U 

where y = {adn^^^),a = 32/(3^/7?), 7 = 27ra/[(67r2)2/3^]. Por any value of u the bulk 
chemical potentials of (E]) and ([6]) are continuous at unitarity. The only free parameter 
u in this expression should be adjusted to have good overall agreement to observables in 
the BEG side of crossover. A set of equations, similar to ([T]) and (E]), for fundamental 
bosons (and not for composite dimers) produced results for energy [15] of a trapped 
condensate in agreement with MG calculations |44] . 

The form ([6]) for the bulk chemical potential has several desired properties. For any 
z/, this chemical potential for the Fermi SF of dimers in the BEG-unitarity crossover, 
by construction, has the following two leading terms of a dilute uniform Bose gas as 
obtained from the model-independent expression for energy per particle derived by Lee, 
Huang, and Yang [251 EE] 

/i(n, a) = {ATTfi^adu/m) 1 + a{n^/^adf/'^ + ... . (7) 

Higher-order terms of expansion ([7]) has also been considered [15] ; the lowest order term 
was derived by Lenz [16]. 

Although ([2]) and ([6|) refer to Gooper-pair induced BGS SF and the dimer BEG SF, 
respectively, these bulk chemical potentials should be interpreted differently. For the 
BGS SF the bulk chemical potential originates from the kinetic energy of Fermi atoms 
put in different quantum orbitals consistent with the Pauli principle discounted for by 
the negative attractive energy due to atomic interaction. For the dimer SF the bulk 
chemical potential originates solely from the repulsive interaction energy among dimers. 

Gonsidering only the lowest-order term of expansion (iTl) in ([1]), appropriate in the 
BEG limit as — )■ 0+, the dimers obey the usual GP equation ^ with the nonlinear bulk 
chemical potential fj.{n,a) = Anh^adn/m. Gonsidering the first two terms of expansion 
([7]), in the BEG limit, one obtains a modified GP equation for dimers written by 
Fabrocini and Polls [17]. But as (and a) increases and diverges at unitarity, the 
nonlinear bulk chemical potential should saturate to the finite universal value of ([3]) and 
should not diverge like the nonlinear terms of the GP equation and of the Fabrocini- Polls 
equation given by ([7]). The chemical potential and energy should not diverge at unitarity, 
as the interaction potential remains finite in this limit, although the scattering length 
a diverges. In the weak-coupling GP limit, the scattering length serves as a faithful 
measure of interaction. But as the scattering length increases, it ceases to be a measure 
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of interaction. Similarly, the bulk chemical potential of ([2]) should saturate, and not 
diverge, at unitarity. 

It has been shown that the DF GP equation ([T]) is Galilean-invariant and equivalent 
to the quantum hydrodynamic equations for dimers [T5l [32] 
dn „ , ^ 

V ■ (nv) = , (8) 



dt 
dv _ 



2m y/n 



H \-U + fi{n, a) 



0, (9) 



if in (fTl) we take 



^(r,t) = V'n(r,t) e''^^'^'\ (10) 
V =hV9/m, (11) 

where v is velocity of hydrodynamic flow of Fermi SF pairs or dimers of mass m and 
6 is the phase of the wave function. The quantum hydrodynamical equations ([8]) and 
(|9]) are new and the usual classical hydrodynamical equations [5] for Fermi SF can be 
obtained by setting ^ = in ([9]). 

In brief, the desired qualities of the present formulation are the following: (i) 
Galilean invariance and complete equivalence to hydrodynamical equations, (ii) correct 
BCS and BEG limits given by (jl]) and ([7]), respectively. There have been previous 
attempts [30l [SI [32l [33] to formulate dynamical equations for a trapped Fermi SF 
valid in the BGS-BEG crossover. The formulations of [301 [3T1 132] had a different mass 
in the phase-velocity relation f[TT]) and hence led to an inappropriate gradient term 
in ([1]) violating Galilean invariance and destroying the equivalence with the quantum 
hydrodynamical equations. The formulations of [301 EI] used mathematical fitting 
functions for the bulk chemical potential which did not satisfy the weak-coupling limits 
dH) and dZD. 



3. Numerical results 



To see how DF equation ([T]) with ([2]) and (|6]) work in practice, we apply it for the 
numerical study of a Fermi SF in an axially symmetric trap. In the first part of this 
investigation we study the density profile and root mean square (rms) size of the Fermi 
^Li SF for the trap parameters and compare our numerical results with experiment 
[34] . In the second part we calculate the frequencies of radial and axial oscillations of a 
cigar-shaped [21] Fermi SF and compare the results with experiments [351 ESI [37] . 

We solve the partial differential equation ([T]) in the BGS-BEG crossover by 
discretizing with the semi-implicit Grank-Nicolson algorithm and employing imaginary 
time propagation for the study of stationary ground states. For the study of dynamics 
we use real time propagation. To this end we use the FORTRAN programs provided 
in [IH]. For discretization we use the space step 0.04 in both radial and axial directions 
and a time step of 0.0004 (both steps expressed in harmonic oscillator units in axial 
direction) . 
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Figure 1. (Color online) The rms size of the Fermi SF for different magnetic field B 
ioi N — Ax 10^ ^Li atoms with the axially symmetric trap described in [34]. Numerical 
results for ^ ~ 0.44, v ~ \; and ^ = 0.35, v =1 are exhibited by lines and experimental 
results are shown with error bars. The vertical line aX B = 0.834 kG shows the position 
of the Feshbach resonance. 



3.1. Stationary Fermi SF 

The experiment of [31] was performed with ^Li atoms in the BCS-BEC crossover using 
a Feshbach resonance at an external magnetic field of 834 G. The results reported 
there are with a radial frequency of = 27r x 640 Hz and an axial frequency of 
{600B/kG + 32)1/2 g^, where B is the external magnetic field in kG. The 
results of the experiment are classified by the magnetic field B, which has been varied 
near the Feshbach resonance at 834.15 G to control the scattering length and achieve 
the BCS-BEC crossover. The parametrization of the external magnetic field B (G) in 
terms of the Fermi-Fermi scattering length a [19] 

300 



-1405ao 



B - 834.15, 



B - 834.15 

1 H 

2500 



(12) 



with ao the Bohr radius, has been used to calculate the scattering length. The results 
reported are with 4 x 10^ atoms or with = 2 x 10^ molecules. 

The DF equations ([T|) and ([2]) valid in the BCS-unitarity crossover produced 
excellent results [TH [19] for energy of a Fermi SF in a spherically symmetric trap, 
when compared to Monte Carlo calculations [2D], at unitarity for 4 to 30 atoms and 
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Figure 2. (Color online) Linear density profile of the Fermi SF for magnetic field B 
= (a) 764 G, (b) 850 G, (c) 809 G, and (d) 882 G, and for iV 4 x 10^ ^Li atoms, as 
calculated by using the DF equation for different ^ and from experiment |34| . 

along the crossover for 4 and 8 atoms for ^ = 0.44. In the BEC-unitarity crossover the 
relevant equations are ([T]) and ([6]) containing the free parameter u. 

Now we study the rms size and density profile of the Fermi SF across the BCS-BEC 
crossover for different values of the magnetic field B, consequently for different values 
of scattering length. After some experimentation to fit the observed density profile and 
sizes [3lj in the BEC side, we find that u = 1 in provides good overall fit and we use 
this value in all calculations here. At unitarity the Thomas- Fermi results for the rms size 
and density profile can be calculated analytically and the best fit to the experimental 
density profile is obtained for ^ = 0.27 [5], [M]. This suggests the use of a smaller value 
of C, (than ^ = 0.44) in the present analysis. Using the present DF equations we also 
find that C, ~ 0.35, rather than C, = 0.44, gives the best overall agreement with the 
experimental density profile for different B (see Figs. [2] below). (The very small value 
for the Bertsch parameter 0.27 obtained in [M] is probably due to systematic error in the 
experimental analysis using the inaccurate position of the Feshbach resonance at 850 G. 
The later analysis used the more precise location of the Feshbach resonance at 834 G, 
giving a more reliable value of the Bertsch parameter [39].) In figured] we plot the rms 
size vs. magnetic field B for ^ = 0.44 and 0.35 together with the experimental results. 
From the general trend of the results exhibited in figure [H it seems that ^ = 0.35 gives 
good agreement with experiment. 

Next we present the results of linear density in Figs. [2] (a), (b), (c), and (d) for 
magnetic field 5 = 764 G (estimated a = 4.474x10^00), 850 G (a = -2.817 x lO^ao), 809 
G (estimated a = 1.52 x IO^Oq), and 882 G (estimated a = —1.041 x lO^ao), respectively. 
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as calculated from the DF equation for different ^ and compare them with experiment. 
(The present estimates for scattering length are slightly different from those presented 
in [33].) From figure |2] (a) we see that at -B = 764 G the linear density calculated with 
^ = 0.44 and 0.35 are practically the same. This is because this value of B is deep into 
the BEG region and away from unitarity and the result is not very sensitive to the value 
of ^. The GP equation for dimers yields good result there [31]. In figure |2] (b) we find 
that near unitarity {B = 850 G) the experimental result prefers a smaller value of ^, e.g. 
^ = 0.3 over ^ = 0.44 and 0.35 in agreement with [5l|3l|. From Figs. |2](c) and (d) we find 
that for B = 809 G and 882 G, ^ = 0.35 should better explain the experimental trend. 
In view of the results reported in Figs. [D and [2], further experimental and theoretical 
studies are needed to see how universal the constant ^ is and how sensitive it is to the 
applied trap, magnetic field and the number of atoms [39]. The approximate solution of 
the Bogoliubov equation as reported in [33] leads to a slightly inferior agreement with 
these experimental data compared to the present model. 

3.2. Oscillating Fermi SF 

Now we study the breathing oscillation of an elongated Fermi SF in the BGS-BEG 
crossover. The frequencies are obtained numerically from the oscillation [HI [35l |36l [37] 
of a cigar-shaped Fermi SF. In our study we take the ratio of the trap frequencies 
^p/^z = 10 and = 40000 atoms. Before we illustrate our results it is appropriate 
to present the analytical results for the problem suggested by Gozzini and Stringari 
[50] in the case of power-law dependence of the bulk chemical potential fi{n, a) as a 
function of n : fi{n,a) ^ n'^ . The collective radial breathing mode frequency is given 
by Vtp/ujp = ^J^il + 7), and the collective axial breathing mode frequency is given by 

^z/^z = y^(2 + 37)/(l + 7). At unitarity and the BGS limit 7 = 2/3, and consequently, 

VLp/ujp = y^lO/3 ^ 1.82574 and QJu, = y^l2/5 ^ 1.54919; at the BEG hmit 7 = 1, and 

consequently, Vtp/ujp = 2 and flz/^z = y^5/2 ^ 1.58114. 

The radial and axial frequencies of a cigar-shaped Fermi SF are calculated by 
studying the time evolution of a pre-formed condensate in a slightly altered trapping 
condition. The radial and axial sizes execute periodic oscillation of the type reported 
in figure |3] (a), from which the radial and axial frequencies are calculated reasonably 
accurately. In the case of small breathing oscillations, these frequencies are found 
to be independent of the perturbation that initiated the oscillation. The normalized 
frequencies flp/oup and ^z/^z are also independent of Uz, ujp and A^, provided that 
A^ is large and the cigar-shaped condition is satisfied: Up » ujz- The results for 
normalized radial and axial frequencies are plotted vs. l/(Ki?a) in Figs. [3] (b) and (c) 
for ^ = 0.44, where up is the Fermi momentum, and compared with different experiments 
and other calculations. (Unlike in Figs. [1] and [21 we could not find a noticable change 
in the frequecies as ^ is changed from 0.44 to 0.35.) In the axially symmetric trap the 
Fermi energy Sp is defined as Sp = {QNY/^h{ijj'^pijJzY^^ and the Fermi momentum Kp is 
defined by Sp = h'^Kp/m. (Recall that A^ is the number of molecules/pairs of mass m.) 
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Figure 3. (Color online) (a) Typical small oscillation of axial and radial sizes of a 
cigar-shaped condensate after an initial perturbation. Normalized frequencies of (b) 
radial and (c) axial oscillation of a cigar-shaped Fermi SF vs. atan([Ki?a]~^) where 
Ki? is the Fermi momentum in the axially symmetric trap. The arrows denote the 
analytic results of Cozzini and Stringari [50] in the BCS (a — > —0) and BEC (a +0) 
limits. DF: present results with ^ = 0.44,1/ = 1; Monte Carlo: from [5l[35l[51]; Hu: 
microscopic mean- field results [52]; Manini: numerical results from the solution of a DF 
equation (31]; BCS: from |5l|35l|51]; Altmeyer: experiment [35]; Kinast: experiment 
[36j : Bartenstein: experiment [37] . 
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The experimental radial frequencies of figure E] (b) were taken from [55] while the 
experimental axial frequencies of figure |3] (c) were taken from [37]. In figure [3] (b) we 
compare the radial frequencies with the Monte Carlo results of Astrakharchik et al. [51] 
as quoted in [3 [35], the mean-field BCS Bogoliubov-de Gennes model calculation of Hu 
et al. [52], and a numerical solution of a DF equation of Manini et al. [3l]. In figure [3] 
(c) we compare the axial frequencies with the Monte Carlo results of Astrakharchik et 
al. [51] , the mean- field model calculation of Hu et al. [52] , and a numerical solution of 
a DF equation by Manini et al. [31]. The numerical calculation of Manini et al. was 
performed with different cigar-shaped traps and different number of atoms from those 
used in the present study. Also, their gradient term in the DF equation was different 
from ours. The BCS and Monte Carlo calculations for the frequencies are performed 
[5T] by using the BCS [40] and Monte Carlo [27] equations of state in a hydrodynamic 
approach. The input for the calculation of frequencies labeled Hu [52] and BCS [51] are 
both based on BCS mean-field approach and lead to very similar results. Qualitatively, 
all the theoretical results have the same trend. However, the curves labeled Hu and BCS 
do not have the over-shooting while approaching the BEC limit as in the experiment 
and the curves based on the DF approach and Monte Carlo equation of state. In the 
BEC-unitarity crossover the present frequencies are very close to the frequencies based 
on the Monte Carlo equation of state. 

4. Conclusion 

We have suggested a generalized time-dependent Galilean-invariant DF equation for 
a trapped Fermi SF valid in the full BCS-BEC crossover, which is equivalent to a 
quantum hydrodjTiamical equation. The equation of state has the correct model- 
independent limiting behaviors near unitarity and near the BCS [25] and BEC [26] 
limits as established by Lee, Yang, and Huang. While applied to the study of size 
and density profile of a cigar-shaped Fermi SF, the present formulation yields results 
in agreement with an experiment on ^Li atoms [31]. We also applied the present DF 
equation to the study of breathing frequencies of a cigar-shaped trapped Fermi SF and 
the results are in good agreement with experiment [351 [37] . 

We also compare the present results with other theoretical findings [IHl [SH [331 
[5T1 [52] . The size and linear densities exhibited in Figs. [D and [2] are in reasonable 
agreement with those obtained from an approximate solution of the Bogoliubov-de 
Gennes equation [33]. The accurate solution of the Bogoliubov-de Gennes equation 
is numerically far more complicated than the solution of the the present DF equation 
and the present simulation yields slightly better agreement to the experimental data 
[M] than the approximate numerical solution of the Bogoliubov-de Gennes equation in 
[33] . The present frequencies for axial and radial oscillations of a cigar-shaped Fermi 
SF are in excellent agreement with hydrodynamical results [51] based on a Monte Carlo 
equation of state [27] and the precise experiment of Altmeyer et al. [S^ but are in 
disagreement with a mean-field BCS equation of state [IQ] and the experiment of Kinast 
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et al. [3H] except near unitarity. The data of Kinast et al. [5H] away from unitarity 
seems to be inaccurate and the overshooting of the frequencies in Figs. [3] above the 
Cozzini-Stringari hmit [50] in the BEC side of crossover is appropriate. More precise 
experiments can determine which of the theories are more reahstic. The experiments 
are easier in the BEC-unitarity region and in this region there is some disagreement 
between two experiments [351 EH] for the radial frequencies reported in figure [3] (b), and 
the large experimental error in the axial frequencies [3^ does not allow us to choose 
between the mean-field calculations [iQl |52] on the one hand and the present and Monte 
Carlo calculations [51] on the other hand. More accurate experiments are needed to 
resolve the discrepancy. 

Although the present approach yields good results for many experimental 
observations, it has its limitation too. For many dynamical problems, such as in rotation 
and vortex creation, breaking of the weakly bound pairs in the BCS-unitarity crossover 
may play a vital role in the phenomenology of the Fermi super-fluid. At sufficiently 
small negative Kj^a, both the Innsbruck and Duke groups observed a breakdown 
of hydrodynamics in the radial breathing mode, presumably because the excitation 
frequency of the collective mode exceeded the pairing energy in the weakly interacting 
regime [36]. Pair breaking is not permitted in the present DF formulation, which may 
lead to inappropriate result in these problems. No such problem should appear in the 
application of the present formulation to a description of density profile, sizes, and 
frequencies of small oscillation as in this paper. 
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